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Abstract 

We study diffraction corrections to the semiclassical spectral density of 
dispersing (Sinai) billiards. They modify the contributions of periodic orbits 
(PO's), with at least one segment which is almost tangent to the concave 
part of the boundary. Given a wavenumber k, all the PO's with length up 
to the Heisenberg length 0{k) are required for quantization. We show that 
most of the contributions of PO's which are longer than a limit 0(k 2 / 3 ) must 
be corrected for diffraction effects. For orbits which just miss tangency, the 
corrections are of the same magnitude as the semiclassical contributions them- 
selves. Orbits which bounce at extreme forward angles give very small terms 
in the standard semiclassical theory. The diffraction corrections increase their 

amplitude substantially. 
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Dispersing (Sinai) billiards were the first class of billiards for which classical chaoticity 
was rigorously proven [1]. These billiard are defined by boundaries which consist of concave 
circular arcs and possibly neutral (straight) segments. Chaos is generated by the divergence 
of neighboring orbits when they bounce off the concave arcs. Whenever a classical trajectory 
is tangent to an arc, the classical billiard map is discontinuous. This singularity in the clas- 
sical dynamics introduces diffraction effects when the billiards are semiclassically quantized. 
Due to the finite wavelength, the domain which is affected consists of a finite neighborhood 
of the classical tangent trajectories, and is called the penumbra (almost shadow). In the 
present paper we shall analyze the effect of diffraction in the penumbra on the semiclassical 
quantization of dispersing billiards. We shall show that it introduces modifications to the 
standard semiclassical trace formula [2], which are of the same order as the leading terms. 
Moreover, it will be shown that diffraction effects are significant because classical PO's which 
are relevant for semiclassical quantization traverse the penumbra with very high probability. 

Diffraction corrections in the penumbra are to be distinguished from contributions of 
creeping trajectories, which are described by the geometrical theory of diffraction [3]. Creep- 
ing corrections were previously included in the trace formula [4], and were shown to be ex- 
ponentially small. Thus, the present theory applies to orbits which fall outside the domain 
of validity of both the standard semiclassical treatment, and its extension which includes 
creeping orbits. 

For the sake of clarity of presentation, our analysis will concentrate on the simplest 
dispersing billiard, whose boundary T consists of a square S and a concentric inscribed 
circle C. We shall refer to it as the Sinai billiard. We introduce the diffraction effects by 
using a variant of the boundary integral method: The eigenvalues of the billiard are those 
real values of k for which the boundary integral equation 



has a solution. In contrast to the standard boundary integral method, we use the circle 's 
Green function Cr c (r, r'), and the integration is carried out along the boundary excluding the 





2 



circle, i.e. along S. The circle's Green function satisfies (A + k 2 )G c (r, r') = —6(r — r') in the 
exterior of C, with Dirichlet boundary conditions on C and outgoing boundary conditions 
at infinity. 

The contribution of PO's which bounce N times off S to the density of states d(k) } is 
derived from (f ), by considering the iV'th term in the multiple reflection expansion 

Im ~^77T / dsi ■ ■ ■ dsN ^O* 1 ' r2 ) " " " (rjv-i, r JV )|^(r JV , n) . (2) 
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The PO contributions are identified when these integrals are evaluated by the stationary 
phase approximation, with each saddle point corresponding to one PO. For this purpose 
G c (r,r') should be approximated by a sum of terms of the form A e tkL ( r > r )^ where A varies 
smoothly with r and r'. 

The circle's Green function is expressed exactly, in polar coordinates, by an angular 
momentum expansion and by Poisson resummation, as G c (r, r') = J2m=-oo ^c m K r ^ r ')i wnere 

G( m )(r, r') = i |_~ d/ (#f (*r<) + S,(kR)H+(k r< )) H+(kr > )e a * e+ ™ ml . (3) 

In this equation, r> (r<) is the larger (smaller) of r and r', Si(kR) = —Hf(kR)/Hf(kR) is 
the diagonal element of the circle scattering matrix, and the angle difference A# = \6 — 0'\ is 
always taken so that < < tt. We assume that r, r' > R + R(kR)~ 1/3 , i.e. the two points 
are not in the near vicinity of the circle. The term G^\r, r') gives the dominant contribution. 
We shall briefly outline the derivation of the standard semiclassical and creeping results, and 
then, the leading corrections in the penumbra will be derived. 

The domain of r and r' in which the standard semiclassical result holds is called the 
illuminated region (see Fig. 1(a)). In this region, the integral for G^°)(r,r') is evaluated 
using the Debye approximation for the Hankel functions [5] , and the integration is performed 
using the stationary phase approximation. There are two saddle points, which relate to the 
two classical trajectories from r to r': One is direct, and the other reflects once from the 
circle. Using the two contributions in (2), one recovers the standard semiclassical trace 
formula. The Debye approximation for Si(kR) fails if (kR — l r ) ^ (^i?) 1 / 3 , where l r is 



(a) 



illuminated 





FIG. 1. (a) The three regions of r' for a fixed r. The direct and reflected trajectories are shown 
for r' in the illuminated region, and the shortest creeping trajectory is shown in the shadow region, 
(b) The geometrical setup in the penumbra. In the upper part r and r' are in the classically 
illuminated region (b^ > R), and in the lower part they are in the classically shadowed region. 

the angular momentum of the reflected trajectory. This sets the limit of validity of the 
semiclassical approximation, and defines the borderline between the illuminated region and 
the penumbra. 

The term G^°)(r,r') gives the exponentially small contribution of a creeping trajectory 
in the shadow region (see Fig. 1(a)). It is obtained by closing the contour of integration 
in (3) in the complex / plane, and summing over the poles of Si(kR) [3]. Using creeping 
contributions in (2) one recovers the results of [4]. This procedure fails if 7^°) < (&i?) -1 / 3 , 
where 7^) is the creeping angle of the trajectory. The above condition defines the borderline 
between the shadow region and the penumbra. A similar evaluation of the terms G^ m )(r, r') 
with ra/0 results in creeping contributions for any r and r'. 
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In the penumbra, neither of the above approximations hold, and one should use specially 
designed methods for approximating G c (r,r'), such as the one introduced by Nussenzveig 
[6]. The expression for G^°)(r,r') is split into two parts, a direct part and a glancing part, 
which are respectively defined by 

G d (r,r') = - / dlH+(kr)H+(kr')e il * e , (4a) 

o JkR 
j rkR 

G g (r,r') = - / dlS l (kR)H+(kr)H+(kr')e UAe 
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+ - / dl (Si(kR) - l)H+(kr)H+(kr')e UAe . (4b) 

8 JkR 

The limits <TiOO and cr 2 oo are directed infinities in the complex plane, defined in [6], such 
that the full contour from <TiOO to cr 2 oo goes through kR and around all the poles of Si(kR). 

The evaluation of the direct part is similar to that of the direct trajectory contribution in 
the illuminated region. The only difference is that here we take into account the boundary 
of the integral at kR. As a result, the semiclassical contribution of the direct trajectory is 
multiplied by the factor (F(oo) — F(v))/y/2i, where F(x) = C(x) + iS(x) is the Fresnel 
integral function, and 

v=(—Y(R-b d ). (5) 

\ TTZZ 1 



The impact parameter of the direct path is denoted by bd, z = yr 2 — b 2 dl z 1 = yr' 2 — b 2 and 
L = z + z' is the length of the direct trajectory (see Fig. 1(b)). The Fresnel factor equals 
| for exact tangency, and goes to 1 (0) when approaching the borderline of the illuminated 
(shadow) region. 

The main contribution to the glancing part (4b) comes from the vicinity of / = kR. In 
the integrand, Si(kR) is replaced by its transition region approximation for / kR [5]. The 
rest of the integrand is evaluated at / = kR, where the Debye approximation is used for the 
Hankel functions. This way, one obtains 

Gfl(r y ) = _^i^L e .^ + L r , + ^o ))j 

47T Ky/ L r L r i 
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where C ~ 0.996193019928e 87r / 3 was obtained by a numerical integration. The length of the 
line segment from r (r'), to the point where it is tangent to the circle, is L r (L r i), and 7^) 
is the directed creeping angle (see Fig. 1(b)). When bj < R the length L r + L r < + R'y^ is 
equal to the length of a creeping path from r to r'. When bj > R the angle 7^) becomes 
negative. 

The contribution of orbits that traverse the penumbra to the density of states d(k) may 
now be found, using (2) and the appropriate approximation for the Green function in each 
of the segments of the orbit. For isolated PO's that traverse the penumbra once, with a 
direct segment, the standard semiclassical contribution (which is 0(k )) is multiplied by 
the Fresnel factor. Thus, the correction to the trace formula will be of the same order 
as the contribution itself. In the case of a single glancing traversal in the penumbra, the 
contribution is proportional to A; -1 / 6 . Due to the new form of the phase and the prefactor 
of G 3 (r,r'), there is no direct relation between the standard semiclassical result and the 
modified one. The former is very small since at glancing, the classical motion is extremely 
unstable. This is manifested in the standard semiclassical result by the fact that G>(r, r') is 
very small for an almost tangent reflection. Replacing it with G 3 (r,r'), which has a larger 
amplitude, the contribution of the PO is typically enhanced. The extension to PO's with few 
segments that traverse the penumbra is straightforward. Since Cr c (r, r') gives a contribution 
in the classically shadowed region of the penumbra, there will also be significant contributions 
of PO's with classically forbidden segments. 

We shall present now a few examples which illustrate the importance of the penumbra 
corrections. This is well demonstrated in terms of the length spectrum 

D(x) = / dkw(k)e ikx d(k), (7) 
Jo 

where w(k) is a weight function. Every PO contributes to D(x) in a small vicinity of its 
length, and this allows us to isolate the contribution of a single PO (or of few orbits with 
similar lengths) and to compare it to the theory. 

We consider the quarter Sinai billiard (see Fig. 2), with side 1 and radius 0.5. The exact 
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FIG. 2. Penumbra corrections of the length spectrum for four selected cases. Solid lines show 
the quantum (exact) length spectrum. All other lines show the deviations from various semiclassical 
approximations (SCA). The orbits considered are shown below each frame. Vertical bars indicate 
locations of unstable PO's, daggers indicate bouncing ball families, (a) Shortest tangent orbit 
(x = 2). Dotted line - SCA including bouncing ball and edge contributions according to [10]. 
Dashed line - glancing contribution also included. Note the logarithmic scale, (b) Double traversal 
of the orbit considered in (a). Dotted line includes 3 penumbra contributions (see text), (c) 
Pair of almost tangent PO's at x ~ 5.10. Dotted line - standard SCA, dashed line - direct term 
included, sparse dotes - glancing contribution also included, (d) Pair of classically forbidden PO's 
at x ~ 5.24, indicated by double dagger. Notation is as in (c). 
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(quantum mechanical) eigenvalues in the range < k < 300 (total of 5667) were obtained 
using the scattering method described in [7]. The classical PO's were calculated using the 
minimum and the unique coding principles [8,9]. To obtain the results for the trace formula, 
we used the four-fold desymmetrization of the circle's Green function. The weight function 
w(k) in (7) was taken as a Gaussian centered around k = 150, and whose width is a = 40. 

We first consider exactly tangent orbits. They always exist at the edge of the one- 
parameter families of neutral PO's ("bouncing balls") that bounce only between the straight 
segments of the billiard. The simplest of these families is of length 2 (see Fig. 2(a)). Accord- 
ing to the analysis in [10], the semiclassical contribution consists of two terms: One that 
describes the contribution of the family (and is 0{k 1 / 2 ))^ and the other that takes into ac- 
count the effect of the edge of the square (and is 0(k )). Comparison with the exact length 
spectrum shows a significant difference. Application of our analysis in the penumbra, shows 
that in addition to the family and edge contributions, there is a term (9(& _1 / 6 ), coming 
from an isolated exactly tangent orbit. This term is obtained from G g (r, r'). The additional 
contribution successfully accounts for the error of the semiclassical approximation. 

The semiclassical result for the double repetition of this family is similar to the single 
repetition, and again shows a significant deviation from the exact result (Fig. 2(b)). Appli- 
cation of our analysis in this case gives three diffraction terms of orders 0(k°), 0(k~ 1 ^ 6 ) and 
(9(& -1 / 3 ), which reduce the error substantially. 

The case of unstable isolated PO's traversing the penumbra, is demonstrated for the two 
PO's of length 5.10 (see Fig. 2(c)). The two orbits are geometrically similar, except that 
one (the direct) passes very close to the circle, and the other (the glancing) reflects in an 
extremely forward direction. The semiclassical amplitude of the direct orbit is multiplied 
by a Fresnel factor, which is 0.71 exp(— 0.23z) for k = k = 150. This accounts for most of 
the semiclassical error. Including the corrected contribution of the glancing orbit, reduces 
the remaining difference by half. 

Next, we illustrate the corrections due to classically forbidden PO's which traverse the 
shaded part of the penumbra (see Fig. 1(b)). The two orbits shown in Fig. 2(d) serve this pur- 
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pose. Their length is 5.24, one of them cuts through the disc and the other creeps around 
it. (If we were to reduce the radius of the disc continuously, these orbits would coalesce 
at R 0.48 and appear as a single tangent orbit). The diffraction contributions of these 
non classical orbits successfully account for the large difference between the exact quantum 
length spectrum and its semiclassical counterpart at the relevant length. Application of the 
geometrical theory of diffraction [4] (including many modes) results in large deviations from 
the quantum mechanical length spectrum. 

The above examples illustrated in detail the effect of diffraction corrections on the con- 
tributions from particular PO's. In the remaining, we provide some arguments showing 
that most of the relevant PO's are affected. One should bear in mind that the borderlines 
of the penumbra are k dependent. Any orbit (which is not exactly tangent), is excluded 
from the penumbra for a sufficiently high wavenumber. Hence, the standard semiclassical 
contribution of any particular PO is valid in the semiclassical limit k — > oo. Moreover, the 
fraction of phase space occupied by the penumbra, is of order (kL)~ 2 / 3 } where L is a typical 
length of the billiard, and this fraction vanishes as k — > oo. Thus one may nai'vly conclude 
that the global effect of diffraction is negligible, which is not true because of the following 
reasons: For an orbit with n segments, define A = minj \lj — kR\ } where lj is the angular 
momentum of the jth segment with respect to the circle's center. The orbit traverses the 
penumbra (at least once) if A ^ (kR) 1 / 3 . Each segment of the orbit has an a-priori prob- 
ability p (kL)~ 2 / 3 to traverse the penumbra. Assuming statistical independence of the 
segments, and homogeneous coverage of phase space by long PO's, the probability that the 
orbit will avoid the penumbra is (I — p) n ~ exp(— n(kL)~ 2 / 3 ). In order to achieve an energy 
resolution of the order of the mean level spacing, one has to include in the semiclassical 
theory all the PO's whose lengths are smaller than the Heisenberg length Lh ~ kL 2 . An 
orbit of this length has njj ~ kL segments. Thus, for all orbits with n^/ 3 < n < njj the 
probability that the orbit does not traverse the penumbra is very small. The exponential 
proliferation of PO's guarantees that most of the relevant orbits are in this range. Thus, 
simple ergodic considerations give the surprising result, that the semiclassical approximation 
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fails for the majority of relevant PO's in the semiclassical limit k — > oo. The contributions 

2 /3 

of orbits whose length is shorter than rijj are described by the standard semiclassical 
theory with a high probability. As k — > oo, the number of these orbits grows, but they 
comprise a smaller fraction of the set of relevant PO's. 

o ^ ■ ^ ■ ^ ■ ^ ■ ^ 1 
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FIG. 3. The coarse-grained distribution of A/(kR), defined in the text. Only angular momenta 
lj of reflections from the circle are considered. The data is based on 19375 PO's which total to 
107029 reflections from the circle. Note the logarithmic scale. 

In Fig. 3 we present the distribution of A for all the PO's in the quarter Sinai billiard 
for lengths between 7 and 10. The distribution is sharply peaked at the minimal value, 
indicating that indeed most PO's include an almost tangent chord. 

The results introduced in this Fetter clearly show that diffraction corrections may play 
a very important role in the semiclassical quantization of dispersing billiards. There are two 
complementary factors which make the diffraction corrections in the penumbra significant: 
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First, the correction of amplitudes of individual orbits is of the same order as the semiclassical 
contribution itself, and hence must be considered prior to higher order corrections (e.g. [11]). 
Second, the corrections should be applied to a large fraction of the relevant PO's. The effect 
of diffraction corrections on an individual energy level, is still to be investigated. 
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